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Abstract 



We present a method to construct infrared-finite amplitudes for gauge the- 
ories with massless fermions. Rather than computing S'-matrix elements 
between usual states of the Fock space we construct order-by-order in per- 
turbation theory dressed states that incorporate all long-range interactions. 
The S'-matrix elements between these states are shown to be free from soft 
and collinear singularities. As an explicit example we consider the process 
e + e~ — > 2 jets at next-to-leading order in the strong coupling. We verify by 
explicit calculation that the amplitudes are infrared finite and recover the 
well-known result for the total cross section e + e~ — > hadrons. 



1 Introduction 



Soft and collinear singularities which arise in gauge theories such as QED and 
QCD are usually dealt with by summing up physically indistinguishable cross 
sections [JJI2]- This results in a cancellation of these singularities for so-called 
infrared-safe observables and therefore reliable theoretical predictions can be 
made for such observables. Even though this is a perfectly valid approach to 
the problem it is useful to investigate the origin of these singularities and to 
explore the possibility of avoiding them altogether. 

The origin of the problem lies in the long-range nature of the interac- 
tions. As a consequence the usual in and out states do not evolve in time 
asymptotically according to the free Hamiltonian. It is this breakdown of 
the standard assumption that results in the non-existence of the scattering 
operator. Thus, if we want to avoid infrared singularities from the outset we 
have to construct true asymptotic states and compute transition amplitudes 
between them j3J|I|. 

This program has been carried out initially for QED with massive fer- 
mions [5| and then many steps have been made to extend it to soft singular- 
ities in non-abelian theories 111 IE] • It is possible to construct asymptotic 
states (generalized coherent states) which include multiple soft gluon emis- 
sion to all orders in the coupling 9J. It can be shown that the S'-matrix 
between such states is free of soft singularities [TUl lllj. Apart from the more 
complicated structure of the soft singularities due to the self-interaction of 
the gauge bosons there is the additional complication of collinear singularities 
in a non-abelian gauge theory. Due to the collinear singularities the asymp- 
totic Hamiltonian is more complicated J21 EH E2] and the prospect of 
being able to include these effects to all orders in perturbation theory are 
not very promising. But the idea of constructing an asymptotic Hamiltonian 
that takes into account the asymptotic dynamics and using the correspond- 
ing evolution operator to dress the usual states [TBI El can still be applied 
and is not tied to any particular theory. In particular, four-point interactions 
that are present in non-abelian gauge theories can be incorporated [TSj . 

In the present work we investigate the practical feasibility of constructing 
asymptotic states whose ^-matrix elements are free from soft and collinear 
singularities. We are not so much interested in general considerations but 
rather try to establish a method to define and explicitly compute infrared 
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finite amplitudes order-by-order in perturbation theory. Apart from the con- 
ceptional advantage of avoiding divergent amplitudes such a method would 
have a variety of practical advantages. Obviously, the finiteness of the ampli- 
tudes would facilitate a completely numerical approach to the calculation of 
amplitudes. This also applies to the combination of fixed order results with 
parton shower Monte Carlo programs. 

After some general remarks about infrared singularities in Section El we 
discuss how to construct the dressed states in Section El This discussion 
will be completely general. In order to test the practicability of our method 
we consider a simple process in Section |3J We compute the infrared finite 
amplitudes for e + e~ — > 2 jets at next-to-leading order in the strong coupling 
and verify that upon integration over the phase space they reproduce the 
well-known total cross section for e + e~ — > hadrons. Technical details of this 
section are relegated to an Appendix. Finally, in Section we present a sum- 
mary and suggestions on how to improve the method to make applications 
to more complicated processes more feasible. 

2 Infrared singularities 

We are concerned with gauge theories with massless fermions. Such theories 
are plagued by infrared singularities. Infrared singularities are related to ei- 
ther arbitrarily soft gauge bosons or arbitrarily collinear gauge bosons and/or 
fermions. 1 More precisely, if we define our external states in the usual way 
by acting with creation operators on the vacuum, then higher-order S'-matrix 
elements between such external states contain infrared divergences. 

Many attempts have been made to define amplitudes that are well defined, 
i.e. do not contain such singularities. However, most attempts were restricted 
to soft singularities. In particular it has been shown that in an abelian 
gauge theory with massive fermions it is possible to define external states 
whose S'-matrix elements are free from infrared singularities to all orders in 
perturbation theory [S]. 

An abelian gauge theory with massive fermions does not contain collinear 
singularities. This simplifies the situation considerably. As soon as we con- 
sider the non-abelian case, however, we cannot avoid the appearance of 

1 We use the term "infrared singularities" for both, soft and collinear singularities 
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collinear divergences. The reason is that in the non-abelian case a mass- 
less gauge boson can split into two arbitrarily collinear gauge bosons. Such 
a splitting results in a collinear singularity. Thus, giving the fermions a mass 
does not protect us from collinear singularities. 

The aim of this work is to investigate the possibility of defining infrared- 
finite amplitudes for a massless, non-abelian gauge theory. The application 
we have in mind is, of course, QCD. In most applications of perturbative 
QCD the quarks (at least the light flavors) are treated as massless. Since we 
will have to deal with collinear singularities anyway in a non-abelian theory, 
we might as well take the common approach and treat the quarks as massless. 

2.1 The conventional approach 

Before we tackle the problem of defining infrared-finite amplitudes in a sit- 
uation where collinear singularities have to be taken into account as well, it 
might be useful to remind ourselves how infrared singularities are dealt with 
in the conventional approach. 

In the conventional approach, which in the following we shall call the 
cross-section approach, we evaluate amplitudes between conventional exter- 
nal states that are obtained simply by acting with creation operators on 
the vacuum. Since such amplitudes contain infrared singularities we regu- 
larize them. In virtually all applications dimensional regularization is used, 
whereby the amplitudes are evaluated in D = 4 — 2e dimensions. The infrared 
singularities then reveal themselves as poles in 1/e. 

In order to relate these amplitudes to measurable quantities, we will have 
to integrate the (squared) amplitudes over the phase space, weighted by a 
measurement function that defines the quantity we are interested in. The 
key point is that we can only obtain theoretical predictions for quantities 
that are infrared safe. An infrared safe quantity is one that does not depend 
on whether or not a parton emits an arbitrarily soft gluon. Also it must not 
depend on whether or not a parton splits into two collinear partons. For such 
quantities the infrared singularities present in the amplitudes are cancelled 
by infrared singularities that appear due to the phase space integration and 
after the cancellation the regulator can be removed. Even though there are 
various general methods available to perform such calculations at next-to- 
leading order JH] the appearance of infrared singularities make the phase- 
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space integration a non-trivial problem. 



2.2 Asymptotic interactions 

One of the basic assumptions in the derivation of Greens functions or tran- 
sition amplitudes is that in the limit t — > ±00 the external states are free. 
Thus, it is assumed that all interactions vanish rapidly enough so that they 
can be neglected in the distant past and future. More precisely, we assume 
that for any exact state vector \^{t)) (in the Schrodinger picture) of the full 
Hamiltonian H we can find a corresponding state vector | of the free 
Hamiltonian H such that 

= e -^|*(0)) -> |$(t)) =e- itHo |$(0)) (1) 

for t — > ±00. We then specify the in- and out-states and ($/| by a com- 
plete set of quantum numbers of H and compute the transition amplitudes 

(S/lfttfi+ISi) = (QflS^i) (2) 

where we expressed the scattering operator S in terms of the usual M0ller 
operators 

ft T = lim e iTH e- irHo (3) 

Since S commutes with H energy is conserved and the amplitude, eq. (J2J), 
is proportional to S(Ei — Ef). This is a point we will come back to below. 

The important point is that in a gauge theory with conventional external 
states the above mentioned assumption is simply incorrect. The interactions 
due to massless gauge bosons are long-range interactions and do not vanish 
rapidly enough for t — > ±00. It is this violation of our basic assumption that 
results in ill-defined, infrared- singular amplitudes. In other words, given an 
exact state \^(t)), in general there does not exist a free state \$(t)) such that 
eq. (JTJ is satisfied. 

If we want to tackle the problem at its root, we will have to compute 
amplitudes between modified external states. The long-range interactions 
that cause the problem will have to be included in the external states them- 
selves. This amounts to replacing our basic assumption with the following: 
we can find an asymptotic Hamiltonian Ha such that for any exact state 
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vector \^(t)) we can find an asymptotic state |H(t)) such that 

\ty(t)} =e- itH \^(0)} -> |3(*)> = e~ itHA \E(0)) (4) 
for t — ► ±00. We then define modified M0ller operators 

Q AT = lim e ^H e -tTH A 

and compute modified S'-matrix elements 

(S/in^JWlSi) = (H/I^IH.) (6) 

Once all true long-range interactions are included in the definition of the 
asymptotic Hamiltonian we are guaranteed that S-matrix elements between 
such asymptotic states are free from infrared singularities |1U1 ITT] . These 
S-matrix elements in turn are related to measurable quantities. In the same 
way as several (squared) amplitudes contribute to a physical cross section in 
the conventional approach, there will be several (squared) S-matrix elements 
between asymptotic states contributing to an observable. The crucial differ- 
ence is that contrary to the cross-section method all these contributions are 
separately finite. 

To find the asymptotic Hamiltonian we have to split the interaction 
Hamiltonian into a "hard" and a "soft" piece 

Hj = H H (A) + H S (A) (7) 

and define i/^(A) = Hq+Hs(A). The separation of the Hamiltonian into two 
pieces is by no means unique. The only requirement is that Hg(A) includes all 
true long-range interactions. Thus, the emission of a soft gauge boson and the 
splitting of a parton into two collinear partons has to be included. But there 
is a lot of freedom on how precisely we make the split between soft and hard 
interactions. To indicate this arbitrariness we include the parameter A in the 
notation. Later, we will often use the abbreviated notation = Hs{A). 
Let us stress that even though we will call = ifs(A) the soft Hamiltonian, 
it includes all long-range interactions that potentially give rise to infrared 
singularities. In particular, the soft Hamiltonian also includes the splitting 
of a parton into two collinear partons. 

In what follows we will define Ha and the asymptotic states more care- 
fully. We then show how infrared-finite S'-matrix elements are related to 
conventional amplitudes. Finally, we consider a simple explicit example, 
e + e~ — > 2 jets at next-to-leading order. 
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3 Infrared-finite amplitudes 



In order to obtain infrared-finite amplitudes we have to find true asymptotic 
states and evaluate the (modified) S'-matrix elements between such states, as 
given in eq. (JBjl. Measurable cross sections then are constructed out of these 
infrared-finite matrix elements. 

Since Sa commutes with Ha we conclude from eq. © that energy is 
also conserved for the modified transition amplitudes. However, given an 
asymptotic Hamiltonian it is generally not possible to find the corresponding 
eigenstates |Hj). First of all, being an eigenstate of the asymptotic Hamil- 
tonian which includes all soft interactions, the true states |Hj) would corre- 
spond to quasi hadrons. In our strict perturbative approach we will never 
be able to describe bound states and, at each order in perturbation theory, 
the states |Hj) will correspond to some sort of "jet-like" states. In particular, 
these states are colored and will have to be related to hadronic states using a 
hadronization model. This is as in the cross-section approach and is an issue 
that we do not address in this paper. Still, we have to relate the asymptotic 
states to conventional free states |<&$) order by order in perturbation 
theory. We are then led to compute matrix elements 

M fi = ($f\S A \$i) (8) 

order by order in perturbation theory and relate them to physical (infrared 
safe) cross sections. It has been argued previously that matrix elements as 
given in eq. (jSJ) are also free of soft singularities (TU1 E]. In the present 
article we extend this to include collinear singularities (see also [E]). For 
more details on this issue we refer to Section 13.51 

3.1 Notation and conventions 

Before we proceed let us fix our notation and conventions. Whereas part 
of the discussion so far was done in the Schrodinger picture we will now 
turn to the interaction picture. Thus all operators and states are now to be 
understood to be given in the interaction picture. 

To start with we construct the usual states of the Fock space 

• ■■'iM'i) ■ --tow ■■■)= n &t fe) n^(Pi) im^) io> (9) 

i j k 
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where b\ (T and denote the creation operators for fermions, antifermions 
and gauge bosons respectively and we suppressed the helicity labels. We will 
generically denote such states by \i) and (/|. Of course, we have to keep in 
mind that the states as given in eq. Q are not normalizable and we tacitly 
assume they have been smeared with test functions. Thus, we are really 
concerned with wave packets. However, we assume they are sharply peaked 
around a certain value of the momentum such as to represent a particle beam 
with (nearly) uniform, sharp momentum. 

The creation and annihilation operators satisfy the usual (anti) commu- 
tation relations 

[a{X 1 ,k 1 ),a\X 2 ,k 2 )} = -(27r) 3 2w(fc 1 )( 7AlA2 8{k 1 - k 2 ) (10) 
{&(ri,£i),&t(r 2 ,£ 2 )} = (2n) 3 2iu(k 1 )5 rir2 5(k 1 -k 2 ) (11) 
{diruhlrffah)} = (2n) 3 2cu(k 1 )5 rir2 5(k 1 -k 2 ) (12) 

with uj(kj) = \ki\. Note that the ordering used in eq. (JHj) implies a certain 
phase convention. Of course, all amplitudes are only defined up to such a 
convention. 

The field operators are given by 

^ Q = Jdk (u a (r, k) b{r, k) e~ ikx + v a (r, k) d+(r, k) e +ikx ) (13) 
V a = Jdk (u a {r, k) tf(r, k) e +ikx + v a (r, k) d{r, k) e~ lkx ) (14) 
= Jdk (e M (A, k) a{\, k) e~ ikx + e*(A, k) a f (A, k) e +ikx ) (15) 

where we defined 

dk = °L (16) 

and the sum is over the two helicities of the fermions or gauge bosons respec- 
tively. 

Once the interaction Hamiltonian is given we can compute the evolution 
operator and obtain in the interaction picture 

U{t, t ) = T exp (-i J* dt #/(*)) (17) 
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The M0ller operators are given by Q± = U(0, =Foo) and, thus, the scattering 
operator S is related to the evolution operator 

s = = u(+oo, o)u(o, -oo) (is) 

This allows us to find the S'-matrix elements between some initial and final 

state 

(f\S\i) = (f\T exp (-i J^dtHftj) \i) (19) 

where \i) and (/| are states as defined in eq. 0. Inserting the explicit form of 
H I into eq. (jTUjl allows us to compute S'-matrix elements. Of course, in prac- 
tice such a calculation is nothing but the computation of the corresponding 
Feynman diagrams. 



3.2 Definition of infrared-finite amplitudes 

In analogy to eq. (|T7|) we define a soft evolution operator 

f/ A (Mo) = T exp (-iJ t <ttH A (tfj (20) 

where we only include the soft Hamiltonian H A (t) = Hs(A,t). Acting on 
a certain state, the soft evolution operator modifies this state by allowing 
for soft and collinear emissions. Then, the usual Feynman-Dyson scattering 
matrix S can be decomposed as 

S = U(+oo, 0)17(0, -oo) = ni_S A (A)Vl A+ (21) 

where we have introduced the soft M0ller operators Qa± = Ua{0, Too). More 
explicitly, we have 

fi A _ = T exp (-i J™ dtH A (tf) (22) 

poo roo rt 

= l-i dtH A (t) + (-i? dt dt' H A (t)H A (t') + . . . 
Jo Jo Jo 

= l-i dtH A (t) + ^—^ dt dt'T{H A {t)H A {t')} + ... 
Jo 21 Jo Jo 

Eq. (}2*Tj) defines a modified scattering operator S^(A). This operator has the 
crucial property that it includes at least one hard interaction and, therefore, 
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matrix elements (f\S^(A)\i) of this operator with ordinary external initial 
and final states as defined in eq. (JOJ) have no infrared singularities. If we 
define dressed initial and final states, and ({/}| according to 

I0» = (23) 

({/}| = (/I«a- (24) 

then 

({/}|5|{<}> = </|fiU(A)|i> (25) 

Thus, the S'-matrix elements of dressed states are free of infrared singularities. 
We should stress that dressed states are not asymptotic states, i.e. they are 
not eigenstates of the asymptotic Hamiltonian. 

Let us look at a dressed final state somewhat more carefully. We obtain 
a dressed final state by acting with Qa- on a final state as defined in eq. • 
We denote this dressed state by adding curly brackets. 

f ({q(pi) . . . q(pj) . . . g(p k ) •••}! = (q(Pi) ■ ■ ■ q(Pj) ■ ■ ■ 9(Pk) ■ ■ ■ \^a- (26) 

Once the asymptotic Hamiltonian is fixed eq. (|26|) is a unique relation, order 
by order in perturbation theory, between an ordinary final state (/| and the 
corresponding dressed final state ({/}|. A similar relation holds for dressed 
initial states. 

\{q{pi) ■ --qijpj) ■ --gipk) ■ ■ •})< = *4+l?(pi) • ■ -?(Pi) • ■■g(pk) ■ ■ •) (27) 

In what follows we will suppress the labels / and % but keep in mind that the 
states \ {q(pi) . --q(pj) ■ ..g(pk) ■ ■ •}> and ({q(pi) . --q(pj) ■ --gijPk) • • -}| are not 
conjugates of each other. Also, we would like to stress that all these states 
are states in the usual Fock space. Of course, this implicitly assumes that we 
use some kind of regularization for the infrared singularities in intermediate 
steps. 

The soft M0ller operators dress the usual non-interacting external states 
with a cloud of soft and collinear partons. Since the infrared behavior of 
i?A and the full interaction Hamiltonian are the same by construction, this 
dressing generates infrared singularities that cancel those generated by the 
full scattering operator. 

There are two main differences between the soft (/collinear) M0ller op- 
erator, eq. and the usual scattering operator, eq. (|19p. Firstly, £Ia± 
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involve only the soft part H A of the interaction Hamiltonian. Secondly, the 
time integration in the soft M0ller operator runs only from to oo rather 
than from — oo to oo. 

The fact that the time integration is restricted to t > is related to 
the loss of Lorentz invariance in the amplitudes A4fi, eq. (jHJ). This is to be 
expected since Sa does not commute with H Q and, therefore, M. fi is generally 
not proportional to an energy conserving S(Ei-Ef). Instead, individual parts 
of the amplitude will have <5-functions with different energy arguments (see 
eq. (|36Jl ). The difference between these arguments determines the amount 
by which energy conservation can be violated in Ai fi and is related to the 
parameter A. In the limit A — > the amount by which energy can be violated 
tends to 0. Thus, the parameter A determines how much the initial wave 
packets are distorted through the evolution with the soft M0ller operators. 
We will come back to these issues in Section 13.71 

3.3 Factorization of modified ^-matrix elements 

We now turn to the question on how to compute the infrared-finite amplitudes 
defined in eq. (I25|) and how they are related to ordinary amplitudes. 

A possible approach is to start from the right hand side of eq. (J33J). 
This would involve using the explicit form of Sa, given below in eq. (|36|) to 
compute the amplitudes. As we argue in Section 13.51 the structure of Sa is 
such that no infrared singularities occur. This opens up the possibility of 
evaluating the amplitudes numerically. We have to keep in mind, however, 
that there are still ultraviolet singularities which will have to be removed 
by renormalization. In order to take an entirely numerical approach the 
renormalization procedure would have to be done at the integrand level (201 ■ 

We will take a somewhat different approach in that we start from the left 
hand side of eq. (|25|) . We relate the infrared finite amplitude to ordinary 
amplitudes by inserting a complete set of states twice 

({f}\s\{i}) = (mA-sni + \i) = <j\n A -\f) ® (f\s\i') ® (28) 

Note that in the final expression all states are ordinary Fock space states as 
defined in eq. In this way, the infrared finite amplitude is split into three 
pieces. First, there is an ordinary S'-matrix element, (f'\S\i'). The other two 
factors are dressing factors for the initial and final state. All these pieces are 
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infrared divergent and only the complete amplitude is infrared finite, order 
by order in perturbation theory. The ultraviolet singularities appear only 
in (f'\S\i') and are dealt with as usual by renormalization. The symbol (g> 
denotes integration over all momenta and summation over all helicities of the 
state under consideration. Thus, for say (f'\ = (q(pi,ri)q(p2,r 2 )g(p3,r 3 )\ = 
{q P iq P 2g P 3\ we have 

\f) ® (f'\ = dpidp 2 dp 3 \q p iq P 2g P z) {q P iq P 29 P z\ (29) 

ri r 2 r 3 J 

We should stress that eq. (|28|) implies that the dressing is not done for each 
external parton separately. The dressing factors (f\QA-\f) do contain terms 
that factorize into separate contributions for each parton, but they also con- 
tain color correlated contributions. 



3.4 Dressing factors 

As we have seen in eq. (128(1 infrared-finite amplitudes are composed of three 
factors. First, there is an ordinary amplitude, (f'\S\i'), computed in the 
usual way using ordinary Feynman rules. Then there are the two dressing 
factors, one for the initial and one for the final state. The calculation of 
these dressing factors is somewhat different from the calculation of ordinary 
amplitudes and it is useful to look at this in some more detail. 

For concreteness we consider the calculation of a final state dressing fac- 
tor. The starting point is eq. ()22j) . Let us stress again that since the time 
integration in eq. (|25|) is from to oo we break Lorentz invariance right from 
the beginning. Of course, in the final result for a physical quantity Lorentz 
invariance will be restored. In fact, the calculation has many features of (old- 
fashioned) time-ordered perturbation theory. Most notably, all particles will 
be on-shell. Three-momentum will be conserved in all vertices, but energy 
will not be conserved. 

A typical term of the (asymptotic) Hamiltonian that gives rise to an 
rz-point interaction has the form 

J dx J f[dki V(h)Q(A) e i*-E"i*i e -i*£"M*0 (30) 

i=l 

where u(kj) = \kj\ denotes the energy of the particles and the sign Oi is posi- 
tive (negative) for incoming (outgoing) particles. V(ki) is made up of creation 
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and annihilation operators, eventually accompanied by spinors and/or polar- 
ization vectors and a certain power of the coupling constant. The range of 
integration over the momenta is restricted to the singular regions. This is 
indicated in the notation by 0(A). The precise form of this function is not 
important at the moment. After performing the dx integration we obtain the 
momentum conserving delta function (2tt) d ~ 1 5( d ~ 1 ^ (J2 &iki). However, since 
the t integration is restricted to t > we do not obtain an energy conserving 
S function. Rather we have to introduce the usual adiabatic factor + > 
and use 

roo j 



dte~ lu)t -> / dte- lu)t e~ m = (31) 

Jo u-m K ' 

Of course, if the t integration was restricted to t < we would have 
o r o , 4 

iujt / i+ — iwt +tO + 



dte-' luJt ^ dte~ tuJt e +m = (32) 

and the sum of eq. (JBTj) and eq. indeed results in 2ttS(uj). 

To summarize, for an n-point vertex in the calculation of a dressing factor 
for a final state we have to use 

/ ndfe,(2 7 r) D - 1 ^- 1 )(E^) - 9 1 A) V&) (33) 

i=i E^(Ki) - «o+ 

Were it not for the 0(A) function and the restriction of the t- integration to 
t > this would lead to the standard Feynman rule. 



3.5 Finiteness of modified S- matrix elements 

In this subsection we substantiate our claim that matrix elements as defined 
in eq. (JSJ) or eq. (|23|) are free from collinear and soft singularities. We start 
from the definition 

s A (A) = n A _sni + (34) 

and use the explicit form of the soft M0ller operator and S to express Sa(A) 
in terms of Ha and iff/. Furthermore, we observe that according to eq. (|30j) 
the time dependence of the Hamiltonian H(tj) is given by 

H H (tj) = hj e -i< * m] ; H A (tj) = sj m ? (35) 
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where Wj = J^^i^ih) is the sum of the energies of the particles associated 
with the corresponding n-point vertex and hj and Sj are time independent. 
Performing the algebra and the t-integrations we obtain up to third order 



Sa 



= 1 — 2m h\ 8 {mi) 



(36) 




h 2m [hi, s 2 ] 

8{mi + m 2 + m 3 ) 



8{mi) - 8{mi + m 2 ) 



+ 2m hi{h 2 + s 2 )h 3 



mim 3 



2m [[hi,s 2 },s 3 ] 




+ 



m 2 {m 2 + m 3 ) 



8{mi) 



) 



— 2m Sih 2 h 3 



8{m 2 + m 3 ) — 8{m\ + m 2 + ro 3 ) 



— 2mhih 2 s 3 



mim 3 

8{mi + m 2 ) — 8{mi + m 2 + cc^) 



First of all we notice that all the purely soft terms sis 2 . . . vanish. This holds 
to all orders and is crucial to ensure that Sa is free from infrared singularities. 
Infrared singularities potentially arise if m^ — > 0. This corresponds to either a 
soft or collinear emission at the corresponding vertex. Let us now go through 
the terms in eq. (j3T?j) and check that for none of them such a singularity can 
occur. For this to be true we have to define hi such that it vanishes for 
m.i — > 0. This can be achieved by choosing the 0(A) in eq. ()3())1 accordingly. 

We start by looking at the second order terms, given in the second line of 
eq. (|3fjj) . The only potential singularity in the first term is mi — > 0. This is 
harmless since h\ = in this limit. In the second term we have the potential 
singularity m 2 — ► which is not prevented by s 2 . However, in this limit the 
term is proportional to 8{mi) and the same argument as for the first order 
term applies. 

The arguments for the third order terms, given in the third to sixth line 
of eq. (jHS|) are similar. The only dangerous limits in the third line term 
for example are m\ — > and m 3 —>■ 0. Both of these are prevented by the 
presence of hi and h 3 . Considering the term in the fourth line, we first note 
that 8 {mi) hi = 0. As a result there is no problem with the limit m 2 —>■ and 
m 2 — > —m 3 . Furthermore, the singularity in the limit m% — > is prevented by 
the presence of hi and the limit m 3 —>■ is made harmless by the combination 
of 8 functions. Similarly, the terms in the fifth and sixth line are finite in the 
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limit vj\ — > and to 3 — > 0. Thus we see that (up to this order) there are no 
singularities in Sa as long as hi is chosen to vanish for w — > 0. 

We mention again that Sa does not only contain terms proportional to 
S(zui + zu 2 + tu 3 ) but also terms with "incomplete" (^-functions. These are the 
energy violating terms mentioned above. We also remark that the absence of 
terms containing 1/ '(wi+Wj) in Sa (the corresponding term in the fourth line 
of eq. vanishes) justifies our initial claim that all infrared singularities 
are related to limits Wi — > 0. 

3.6 Construction of infrared finite amplitudes 

The expression given in eq. (J28|) is a (double) sum over all possible interme- 
diate states \f')(f'\ and However, if we compute the amplitude to a 
certain order in the coupling constant, only a very limited number of inter- 
mediate states contribute. It is for example clear that at order O(g ) the 
dressing factor (f\D,^\f) is zero, unless / = /'. From this we see that at 
leading order in perturbation theory the amplitude is the same 
as (f\S\i). 

Including higher-order corrections this identity will, of course, not hold 
any longer. At order 0(g l ) the states / and /' can be different. To get a 
non-vanishing contribution they must be related either by adding a (soft or 
collinear) gluon or by exchanging a quark-antiquark pair by a gluon. 

In order to illustrate this in more detail, let us consider a concrete pro- 
cess. To simplify matters we consider with no partons in the initial 
state. What we have in mind is for example the process e + e~ — > 7 — ► jets. 
As long as we treat this process at leading order in the electromagnetic cou- 
pling but at higher order in the strong coupling, g, we encounter only final 
state singularities. Thus, for the purpose of understanding how the dressing 
removes the infrared singularities we can restrict ourselves to the final state 
partons and treat the initial state simply as |0). 

Before writing down eq. (|28j) more explicitly for the process under consid- 
eration, let us introduce a somewhat more compact notation. We will denote 
the momenta and helicities of the partons in the intermediate state /' by qi 
and Si respectively and use the notation q qi = q(q iy Sj) etc. The momenta 
and helicities of the partons in the final state / on the other hand will be 
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denoted by Pi and r$ and we use q pi = q(Pi,ri). The order 0(g n ) terms of 
the dressing factors are then denoted by 



g n W (n) (g p i, q p2 , g p3 ---; g q i, q q i ■ ■ •) 



(37) 



(q(Pun)q(p2,r 2 )g(p3,r 3 ) . . . \tt A _\q(q u Si)?(g 2 , s 2 ) • • •} „ 

9 

Similarly, we denote the order 0(g n ) terms of the amplitude by 

g n A n \q{q u Si), q(q 2 , s 2 ),g(q 3 , s 3 ) • • • ; 7) = (38) 
(q(qx, Si)g(g 2 , s 2 )g(q 3} s 3 ) • • • l 5 ! )^, = #™ -4 (n) (?gi; 9g2, <7«3 ■ ■ ■ 5 7) 

and we introduce a notation for the infrared finite amplitudes 

g n A^ n \{q 1 (p 1 , n), g 2 (p 2 , r 2 ), ^ 3 (P3, r 3 ) . . .}; 7) = (39) 
({g(pi,ri)g(p 2 ,r 2 )^(p 3 ,r 3 )...}|5|0) n = <^)({ %1 

> 9p2, 9p3 ■ ■ ■ };t) 

We always make use of the convention that the helicity associated with mo- 
mentum ^ is 7~j whereas the helicity associated with momentum is Sj. 

Let us now use eq. (|28j) to write down the infrared finite amplitude 
({q(pi,ri)q(p 2 ,r 2 )}\S\0) order by order in perturbation theory At leading 
order we have 

^°\{q pl} q p2 };i) = ({ g (p 1 ,r 1 )g(p 2 ,r 2 )}|5|0)| g0 (40) 

= Vy (0) (g p i, q p2 \ q ql , q q2 ) ® ^ (0) (9 9 i, 9g2; 7) 
= ^ (0) (9pi,9 P 2;7) 

where in the last step we used 

W (0) (g P i,g P 2;g 9 i,g g2 ) = (41) 

(27r) 3 2^(p 1 )5 riSl (5(pi - q x ) {27r) 3 2uj{p 2 )5 r2S2 5{p 2 - q 2 ) 

Eq. (jH)) is simply obtained by noting that f^- = 1 at O(g ), eq. (J22J, and 
using the ( ant i) commutation relations eqs.(fTU|). (fTT|) and (JTJJ). 

At 0(g) the amplitude is zero because for every intermediate state /' 
either the dressing factor (/|^a-|/') or the amplitude (f'\S\0) vanishes. 
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(42) 



= W (0) (g P i, %2\ q q i, q q 2) ® A {2 \q ql , q q2 ; 7) 

+ W (2) (<? P i, q p2 ; q g i, q q 2) ® A (0 \q ql , q q2 ; 7) 

+ >V (1) (q p i, q p2 ; q q i, q q2j 0,3) ® A {1 \q ql , q q2 , g q3 ;>y) 



Pi n 




P2 r 2 



Figure 1: Cut diagrams for 2-particle intermediate state. The term 
W(°) ® A® of eq. corresponds to j = 0, i = 2 and W (2) (8) -4 (0) 
corresponds to j = 2, i = 0. 

The first term on the right hand side of eq. (}4*2"j) is nothing but the 
usual one- loop amplitude multiplied by the O(g ) dressing factor and, us- 
ing eq. (|4"Tj). can be written as 



The second term is also a two-particle cut term, but this time it is the usual 
tree-level amplitude multiplied by the next-to-leading order dressing factor. 
These two terms are shown in Figure ^ 

The third term in eq. (J42j) is of a somewhat different nature as it is a 
three-particle cut diagram, as illustrated in Figure 121 The dressing factor 
W (qpii %2\ q q i, q q 2, g q 3) is zero unless the gluon g q3 is either soft or collinear 
to the quark or antiquark. Thus, the dressing factor projects out the infrared 
singular piece of the bremsstrahlung amplitude. This is exactly the piece that 
is needed to render the full amplitude ^4^ 2 ^({g p i, q p2 }; 7) finite. 

In the next section we will calculate this amplitude explicitly and check 
that the infrared singularities present in the three terms of eq. (142)1 cancel in 
the sum. 



>V (0) (g P i, q P 2\ q q i, q q 2) ® A (2) (q ql , q q2 ; 7) = A {2) (q pl , q p2 ; 7) 



(43) 
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pi n 




P2 r 2 



Figure 2: Cut diagram for 3-particle intermediate state. 

The construction of the amplitude at higher orders in g follows the same 
pattern. For any odd power of g the amplitude vanishes for the same reason 
as it vanishes at 0(g). At 0(g 4 ) it is given by 





w (0) 


(ipl 


%2 


Qql 


q q 2) ® A {4) (q qU q q2 ;-f) 


+ 


w (2) 


(q p l 


%2 


Qql 


Q Q 2) ® A {2) (q ql ,q q2 ;7) 


+ 


w< 4 > 


(ipl 


%2 


Qql 


q q 2) ® A (0) (q ql ,q q2 ;i) 


+ 


w « 


{qpl 


%2 


Qql 


q q 2, g q s) ® A {3) (q ql} q q2} g q3 ; 7) 


+ 


w (3) 


(qpl 


%2 


Qql 


q q 2, g q s) ® A {1) (q ql} q q2} g q3 ] 7) 


+ 


w (2) 


(ipl 


Qp2 


Qql 


q q 2, g q 3, g q i) A (2) (q q i, q q2 , g q 3, g q \\ 7) 


+ 




(qpl 


Qp2 


Qql 


q q 2, q Q 3, q q i) ® A {2) (q ql , q q2 , q q3 , q q4 ; 7) 



The separate terms in eq. (|44|) are infrared divergent but in the sum all these 
divergences cancel. This can be seen by looking at a particular Feynman 
diagram, for example the one shown in Figure El and realizing that eq. (|44j) 
is nothing but the sum over all possible cuts. Since the dressing factors 
are constructed such that in the infrared limit they correspond to the usual 
amplitudes it is clear that the infrared singularities in ^4^ 4 ^({g p i, q p2 }; 7) have 
to cancel in the same way as they cancel in ordinary cut diagrams. The 
first term of eq. (|44j) corresponds to the ordinary two-loop amplitude and 
is represented by cut 1. The other two-particle cuts, the second and third 
term, are represented by cut 2 and 3. There are two three-particle cut terms, 
term 4 and 5. Finally, for the diagram under consideration, there is one four- 
particle cut contribution, namely term 6. For a certain Feynman diagram 
not all terms of eq. (J44)) are present. In our case, the last term of eq. (|44j) 



A (4) ({q p i,q p2 }]"f) 



(44) 
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which is another four-particle cut contribution is missing. 




1 6 4 2 5 3 



Figure 3: All possible cuts of a Feynman diagram representing the 
various terms of eq. (|44|) . 

We should stress that our approach to construct infrared finite amplitudes 
is by no means restricted to amplitudes with final state singularities only. 
Initial state singularities are dealt with by dressing the initial state, as can 
be seen in eq. ()28|). 

In fact, the dressing of the initial state would even be needed for processes 
as discussed above, i.e. with say only a 7 in the initial state. Above and in 
the rest of this paper we have excluded any QED vertices from the soft 
Hamiltonian even though there is a potential collinear singularity at this 
vertex. We do this because we treat the incoming photon as off-shell and so 
it will generate no infrared singularities. 

If we were to include such vertices in the soft Hamiltonian then we would 
generate many more diagrams with non-vanishing initial-state dressing fac- 
tors such as W(q(qi, Si), q{q2, S2), g{q\, S3); 7). We would find though, that all 
such extra contributes would cancel as all diagrams with purely soft vertices 
cancel as described in Section 1331 

3.7 From amplitudes to cross sections 

Once the infrared finite amplitudes have been computed, they can be used 
to compute cross sections for observables related to these amplitudes. The 
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procedure to obtain cross sections from amplitudes depends to some extent 
on the external states we use and deserves some further considerations. 

In the cross-section approach we usually deal with amplitudes that are 
proportional to a four- dimensional delta function. Upon taking the absolute 
value squared, this leaves us with the problem of interpreting the square of 
a delta function. Usually this is dealt with in a rather non-rigorous manner 
by putting the system in a four dimensional box of size V ■ T. The square of 
the delta function is then interpreted as V ■ T times a single delta function. 
The factor V ■ T is cancelled by taking into account the normalization of the 
states and the flux factor, which leaves us with a cross section proportional 
to a single four-dimensional delta function, expressing conservation of four 
momentum. 

The appearance of the square of the delta function is of course related to 
the fact that we usually work with non-normalizable states with a sharp value 
of momentum and energy. In a more rigorous treatment within the cross- 
section approach the in and out states would have to be written as wave 
packets, sharply peaked around a certain value of momentum and energy. It 
can then be shown that the spreading of the wave packet during the scattering 
process can be safely neglected [2T]. As mentioned in Section I2~T1 the precise 
definition of the measurable quantity is given in terms of a measurement 
function. This is a function of the partonic momenta. If we are dealing 
with wave packets rather than sharp-momentum states, the measurement 
function has to be defined in terms of these wave packets. However, as long 
as we deal with wave packets whose spread is well below any experimental 
resolution, we can simply use the normal measurement function with the 
partonic momentum replaced by the central value of the wave packet and we 
get the same result as in the above mentioned, less rigorous approach [2*T] . 

Let us now turn to the situation we encounter if we work with infrared- 
finite amlitudes, defined in eq. (jSJ). As mentioned before, the amplitude is 
then not proportional to an energy conserving delta function, even if we 
were to start with the usual non-normalizable states. Following the proper 
treatment with wave packets, we think of the states \i) and (/| (or 
and ($/|) as sharply peaked wave packets. The states \{i}) and ({/}| as 
defined in eqs. (J23|) and (J24j) are also wave packets. Through the action of 
the M0ller operators, their spread is larger than the spread of \i) and (/| 
and depends crucially on the parameter A. If we choose A small enough 
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such that the spread of the wave packets related to the states \{i}) and ({/}| 
is still smaller than any experimental resolution, we can still compute any 
measurable cross section by using the standard measurement function with 
the partonic momenta replaced by the central value of the wave packet. 

The important point is that we must be able to express any measurable 
quantity in terms of the states and ({/}|. However, since the states 
\{i}) and ({/}| differ from \i) and (/| only by soft and collinear interactions, 
this is nothing but the requirement that the quantity we are dealing with 
is infrared safe. Indeed, the requirement of infrared safety states that the 
quantity must not depend on wheter or not a parton emits another arbitrarily 
soft or collinear parton. But in the limit A — > the soft M0ller operators do 
precisely this. Thus, choosing A small enough ensures that the construction 
of the measurable quantity in terms of the partonic momenta is not affected 
by the soft M0ller operators. 

This solves the problem on how to obtain differential cross sections, once 
the infrared-finite amplitudes are known, in principle. In practice, the explicit 
implementation of this programme is far from trivial and requires further 
investigations. We mention for example that choosing A very small might 
result in numerical problems, similar to the so called binning problem in the 
standard approach. If, on the other hand, we choose A too large (relative to 
the experimental resoluion) the infrared-finite amplitudes are too inclusive 
to allow the computation of any possible physical quantity. It has been 
advocated before that the most convenient choice of Ha is the one that 
precisely corresponds to the experimental resolution While this might 
be true in principle, we think that this is not a practicable way to proceed, 
since then the asymptotic Hamiltonian would depend on the details of the 
experiment. 



4 An example e + e — ► 2 jets at NLO 

We consider the process e + e~ — > j(P) — > 2 jets. At leading order there is 
only one partonic process that contributes, e + e~ — > qq. However, at next- 
to-leading order there is also the process e + e~ — > qqg. Since the initial state 
does not interact strongly we can restrict our considerations to the process 
7 *(P)^2jets. 
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4.1 The conventional approach 

In the conventional cross-section approach we compute the amplitudes for 
the two partonic processes 

A{q p ±,q p2 ;i{P)) = 

A^(q pl ,q P 2-,l(P)) + g 2 A^(q pl ,q p2]1 (P)) + 0(g 4 ) (45) 

and 

A(q pl ,q p2 ,g P 3-n(P)) = gA (1 \q pl ,q P 2,g P 3n(P)) + 0(g 3 ) (46) 
Upon squaring the amplitude and integration over the phase space we obtain 



da = da + g 2 da qq + g 2 da qqg + 0(g 4 ) (47) 



where 



da, ~ \A^(q pl ,q p2 -n(P))\ 2 (48) 
da qq ~ 2Re[A i0 \q pl ,q p2 ;i(P))A {1 >(q pl ,q P 2-,l(P))} (49) 
da qqg ~ |^ (1) (9 P i,g P 2,^3;7(^))r (50) 

The virtual cross section, <ia g q and the real cross section, da qqg both contain 
infrared singularities and only when combined to form an infrared safe ob- 
servable do these divergences cancel. For the total cross section for example 
we obtain 

^a 2 cm e 2 q N c 

*o = 3s (51) 

„ a s /2 3 19 tt 2 \ 
a„ = * C F -Cr _ + _ + _-_ (52) 



< a s ( 2 3 7T 2 \ 

= ^oC F — c r l-^- — -4 + — ] (53) 



where c r = 1 + (9(e) and C F = (iV 2 - 1)/(2JV C ) =4/3. Thus, at next to 
leading order the total cross section is given by 

o"i = cr qq + a qqg = a ^1 + ^ 3C F ^) (54) 
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4.2 The infrared finite amplitudes 

In terms of infrared finite amplitudes, at next-to-leading order a cross section 
is also made up of two contributions. The two amplitudes that contribute are 
those with the final states ({q p iq p2 }\ and ({q p xq~ p2 g p z}\- However, the crucial 
point is that both these amplitudes are infrared finite. Up to the order in g 
required they are given by 

A({q p u q p2 }; 7) = ({q P iq P 2}\S\0} = (55) 
J dgidg 2 W (0) (g P i, q p2 , q q i, q q2 ) x A (0) (q ql , q q2 ; l{P)) 

+ J d^dq 2 g 2 W ( - 0) {q pl ,q p2 ;q ql ,q q2 ) X A {2) (q ql , q q2 ;-f(P)) 
+ J d^dq 2 g 2 W (2) (q pl ,q p2 ;q ql ,q q2 ) x A (0) (q ql , q q2 ; 7(P)) 
+ J dqtdq 2 dq 3 g 2 W {1 \q pl ,q p2 ;q ql ,q q2 ,g q3 ) 

x A^uq^g^MP)) + 0(g A ) 

and 

A({q P uq p2 ,g pS };-f) = ({q P iq P 29p3}\s\Q) = (56) 

J dgidg 2 q p2 ,g P 3] q q i, q&) x A {0) (q q i, q q2 ; 7(F)) 

) <?p2, <?gl; <?g2, 5^3 J 

where a sum over the spin/helicities of the intermediate particles is under- 
stood to be included in / d(&, eq. (fTEj). 

4.3 The asymptotic Hamiltonian 

Before we can proceed with the calculation of the infrared finite amplitudes 
we have to define the asymptotic Hamiltonian if a- Once we have if a we 
can obtain the M0ller operator, eq. (J22J), and use it to construct the dressed 
states, eq. (J2EJ), order by order in perturbation theory. 

The only condition on is that it includes all long-range interactions 
from the original Hamiltonian. In order to separate these long range soft and 
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collinear emission terms from the hard emission terms we need to introduce 
(at least) one parameter which we denote by A. The dependence of the 
asymptotic Hamiltonian on this parameter is indicated in the notation 
Once these terms are incorporated into the asymptotic Hamiltonian we are 
free to include any other terms from the original Hamiltonian that we wish, 
as these will only produce finite A-dependent contributions to the two final 
amplitudes. It is clear from eq. (|28j) that for the final result this A dependence 
has to cancel. 

In our case the only term of the interaction Hamiltonian we wish to 
include in Ha is the quark gluon interaction vertex, 

Hj = g J dx : ^T a Y^ ■ A* (57) 

Using eqs. (jl3|14j) and (JT5|l we see that Hj consists of eight terms 

r ~ ~ ~ 8 - - 

Hj = gT a I dhdk 2 dhJ2 V i(^ k ^h) (58) 

i=l 

x exp {-it <7ijv(kj) J <5 (D_1) (^J^Vijkj 
where (suppressing the helicity and color labels) 

V 2 = b\k l )d\k 2 )a{h)-u{k{)^k 3 )v{k 2 ), 

V 4 = -d^k^dik^aik^-^hMhMh), 

V G = dihXhWfy-vikJffyufy, 

V 7 = b\k 1 )S{k 2 )a\k^)-u{k l )f{k z )v{k 2 ), 

V 8 = -d\k 1 )d{k 2 )a\k z )-v{k 2 )f{k^k 1 ) (59) 

The sign factors are +1 (—1) for incoming (outgoing) particles. 

As we only have to include terms in that contribute in the singular 
regions we are free to exclude the Vi in eq. (|59|) for which c'iUj(ki) can never 

— # — * 

equal zero with all u(ki) > such that not all of the u(ki) = 0. We can see 
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from eq. (|33|) that such terms will always be finite. From the remaining terms 
we choose only those that give a singularity in the physically relevant soft or 
collinear regions. This means that for our example we can exclude V 2 and 
Vq, as these only go singular when the two incoming or outgoing quarks from 
the vertex are collinear. We emphasize that for more general processes these 
terms have to be included in the asymptotic Hamiltonian. 

We can confine the remaining terms even further as we are free to choose 
the form of the finite part of Ha- We restrict the integration of the momenta 
fa, fa and fa to just the potentially singular regions. This restriction is 
achieved here by including a theta function, 0(Aj(&i, fa, fa)) in each Vi from 
eq. (jo^jl which will appear in Ha- The form of Aj(A?i, fa, fa) is completely 
arbitrary as long as —>■ 1 in the soft and collinear limits. 

The form of the function that we will take for this example is, 

0(Ai(fci, fa, fa)) = 0(A - I ^Mfa)\)- (60) 

i 

This choice of A(ki, fa, fa) is particularly appropriate because as we see in 
eq. (}3"3"j) . J2j a ij UJ {kj) is the exact form that the singular terms take. This 
theta function therefore restricts the integral to just the regions close to 
these singular limits. 

By splitting up the covariant vertex into pieces and restricting the inte- 
gration to just the singular regions we are removing the manifest Lorentz and 
gauge invariance from the amplitudes. Physical observables will though be 
Lorentz and gauge invariant as we are effectively just performing a unitary 
transformation (as we have regulated the Q± operators) on a known Lorentz 
and gauge invariant result. 

To summarize, for our asymptotic Hamiltonian we take just the vertices 
Vi, I/4, V 5 and V 8 , giving, 



Ha = 9 dkidk 2 dk 3 \ v i(h, fa, h) exp I -it Cijv(kj 

1=1,4,5,8 I \ j=l 

t (D ~ l) (x>y%j e(A- lE^ft)!)}- (6i) 
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4.4 Diagrammatic rules for the asymptotic regions 

We can now take eq. (JoTj) and use it in eq. (J22J) to form the asymptotic 
operator. We could then go on to calculate the dressed states of eq. ()37|) 
with this operator defined between suitable in and out states by using the 
commutation relations between a, b, d, a', b', d> and time-ordered perturba- 
tion theory. However it can be shown that there are a set of diagrammatic 
rules for the asymptotic region which behave in a similar way to Feynman 
diagrams in normal perturbative field theory. Using these we rules we can 
simplify the calculation. 

These diagrammatic rules consist of vertex and propagator 'like' objects, 
but unlike normal Feynman diagrams we must take all time orderings of the 
vertices into account. This is because we base the evaluation of the amplitude 
in the asymptotic region on time ordered perturbation theory. As mentioned 
before energy is not conserved at each vertex and since the range of the time 
integration in the M0ller operators is from to oo there is no overall energy 
conservation. 

As there is a time ordering to the vertices we have both absorption and 
emission rules. These are defined in Figure 0] with time flowing from right to 



We form propagator 'like' objects from the spin sums of fermion spinors 
and an associated energy denominator. Although they are not real propa- 
gators in the normal field theory sense of inverted off-shell two-point Green 
functions, they do represent the transition from one vertex to another. The 
rules for these are shown in Figure El where p M = (p°, —p). 



left. 






Figure 4: The diagrammatic rules for vertices. 
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- mp) \ g " u m ) ' 

diagrammatic rules for propagators. 

As with ordinary field theory we must integrate over all internal momenta 
and so for each propagator in the asymptotic region we must integrate over 
its momentum / d D ~ 1 p/ (27r) D_1 in D — 1 = 3 — 2e dimensions. The rules for 
external particles are exactly the same as for QED or QCD and so do not 
need to be reproduced here. Finally we must include a factor of 1/n! with 
each diagram, where n represents the order in the coupling in the asymptotic 
region. 

As stated before the soft M0ller operators are not necessarily gauge in- 
variant and neither Lorentz invariant. Infrared singularities though will only 
occur in the region where w = J2<JiUj(ki) = 0. In this limit Lorentz invari- 
ance is restored and so the structure of the singularities will also be Lorentz 
invariant. Given that our amplitudes will not be gauge invariant, we will 
perform all calculations including the second term of the gluon propagator 
which ensures that we sum over physical polarizations only. 

4.5 The amplitude A({q(pi), q(p2)}', 7) 

Let us start with the amplitude A({q p i, q P 2}] 7) given in eq. This ampli- 
tude consists of four terms and we will look at each of them in turn. The first 
and second term will be dealt with in Section 14 . 5 . II and 14 . 5 . 21 respectively. For 
the third term of eq. (|55jl we need the dressing factor W^(<7 p i, q P 2, q g i, q q 2)- 
There are various combinations of interaction terms of the asymptotic Hamil- 
tonian that give rise to non- vanishing contributions to W^. These can be 
found using the diagrammatic rules of Section 14.41 We find that there are 
four contributing diagrams. These four diagrams can be split into two classes, 
two self-interaction terms and two one-gluon exchange terms. We consider 
the former in Section 14.5.31 and the latter in Section 14.5.41 Finally, the last 
term of eq. (|53j) will be computed in Section 14.5.51 



V 



a 



Figure 5: The 
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4.5.1 The Born term 



The first term is 

J dq x dq 2 W^(g P i, q p2 ; q qX , q q2 ) x A®(q qX , q q2 ; l{P)) (62) 

and is of order g°. As discussed previously, eq. (I4U|) . this term corresponds 
precisely to the tree-level amplitude. 

A^({q ph q p2 };i(P)) = (^e)5 lJ (p 1 \ 1 a \P2)(27f) D 5^(P-p 1 -p 2 ) 

= A^{q pl ,q p2]1 {P)). 

We use a notation where (pi\ represents the spinor of a massless outgoing 
fermion with momentum pi and similarly \pj) represents the spinor of a mass- 
less incoming fermion of momentum pj. Of course, these spinors depend on 
the helicity of the fermion, but we suppress this dependence in the notation. 
The delta function as usual ensures energy-momentum conservation for the 
process and the 5^ represents the color flow through the diagram. 

4.5.2 The virtual term 

In the same way we see that the second term of eq. (J55j) corresponds to the 
one-loop amplitude 

A^(q pl ,q p2 ; 1 (P)) = (63) 

*(£)(£H44-* + f + s)^***>- 

The infrared singularities appearing in eq. ()63|) will be cancelled by infrared 
singularities of the third and fourth term of eq. (|55|). We should mention 
that the finite term in eq. (}63j) depends on the regularization scheme used. 
The result in conventional dimensional regularization is obtained by setting 
cr = whereas in dimensional reduction we set Cr — 1. 

4.5.3 The self-interaction terms 

The two self-interaction terms are obtained by taking the interacting terms 
Vi, V 5 and tyj, V§ of the asymptotic Hamiltonian as given in eq. (|59p. Since 
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there is a symmetry between these two contributions we only need to calculate 
one of the pair of diagrams. The self interacting term resulting from the 
vertices V\, V5 is shown in Figure El and is given by 



,{2,0} 
l 15 



J d qi dq 2 g 2 MsiQpi, %2\ q 9 i, q&) x A {0 \q ql , q g2 ; l{P)) (64) 
J d ^ d q2 (2n) D - 1 ik kj 



2 J * (2tt 



+ ™7 - V (Pil7^ 2 7^i7 a |P2) 



9(A - |qj(g 2 ) + ^(g 3 ) - ^(pi)|) 6(A - \uj{q 2 ) + a;(g 3 ) - 
2u(q\) (uj(q 2 ) + uj(q 3 ) - a; (pi)) 2a;(g 2 ) (w(<&) + w(<f 3 ) - w(<fO) 
^(fc + 93 - Pi) ^ D ~ l \q 2 + q 3 - qi ) (2n) D 6^ D \P -q x - p 2 ). 

Note that this expression contains a D-dimensional delta function com- 
ing from and two (D — l)-dimensional delta functions coming from 
3-momentum conservation of the vertices in the dressing factor. 




Figure 6: Cut diagram for self interaction with 2-particle intermediate 
state. 



We now proceed to perform the integrals over qi and q 2 , removing the two 
(D — 1) dimensional delta functions. There is an important subtlety here. 
Since the delta functions are (D — 1) dimensional, only the spatial part of 
the 4-vectors is altered. All 4-vectors in the asymptotic region though must 
be on-shell and so we are forced to modify the energy component of these 
4-vectors to preserve this property. Although these modified 4-vectors are 
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4 component objects they no longer transform as tensors. This is simply 
a manifestation of the breaking of Lorentz invariance that occurs in time- 
ordered perturbation theory. To denote such objects we place curly brackets 
of the type { } around them, i.e. we define 



{pi-?3} = (w(^i-S)>R-^) (65) 



We then have 



a 



(2,0, = J_j^ TS j. ;A2w) n sm{p _ pi _ p2) (66) 



/ 



2 

dg 3 [-r + q %2) t3 ) 0(A " ~ m 



2u(p 1 )2uj(p 1 - q 3 )p(q 3 ,Pi - q 3 ) 2 

where we defined 

p(ki, k 2 ) = u}(ki) + uj(k 2 ) - v(ki + h) (67) 

This diagram contains infrared singularities coming from the region where 
q 3 is soft and/or collinear to p\. We discuss its evaluation in the Appendix. 
Multiplying by two to take into account both of the self-interaction diagrams 
we get the final result 

2ag'° } = 2 j ' dq l dq 2 g 2 W^{q pl ,q p2] q ql ,q q2 )yiA^\q ql ,q q2]1 {P)) (68) 



+ J dg 3 6(A- \p(q 3 ,Pi-q 3 )\)fi(pi,p 2 ,q3)jA^(q p i,q p2 ;j(P)) 



with 



7 5 /A\ 2 7tt 2 



7 /A\ 1 /A x 2 

2 + C^-) + 2 



, g (l) 



+ '°S 2 (f)+2'°S 2 (l + |)+4L, 2 (^) (69) 

and where fi(pi,p 2 , q 3 ) is a function that is free from singularities when 
integrated over dq 3 . The explicit form is given in eq. ([99)1 . Note that we took 
care to sum over the physical polarizations of the gluon only and evaluated 
the diagram in the center-of-mass frame p\ = —p 2 . 
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4.5.4 The one-gluon exchange terms 



We now look at the one-gluon exchange diagrams. There are two such di- 
agrams, one for each time ordering of the two vertices. One diagram is 
obtained from taking the vertices Vy, Vg of eq. (JoTJj) . the other from taking 
the vertices V4, V5. These diagrams are symmetric under exchange of all mo- 
menta and so we need only calculate one of them. The diagram shown in 
Figure [7| gives us 



41 ° } = / dgidg 2 g 2 W[f(q pl , q p2 ; q q i, q q2 ) x A {0) (q ql , q q2 ; 7 (P)) (70) 

f lD-l n f lD-l n _^ 93_ T a T b 

a 9l« 92 , _x r>_i 1 ik 1 k 



-ie)g 2 f , D _i d D l q 3 



M9s) V (9s9 3 ) 
9(A - |a>(gi) + o;(g 3 ) - ^(pi)|) 9(A - |cu(p 2 ) + u(q 3 ) - u(g 2 )|) 
2w(?i) + ^(93) - w(pi)) 2u;(g 2 ) (w(&) + w(<f 3 ) - £*>(<&)) 

5 (D - X) (?i + 93 - Pi) S^(p 2 + q 3 - q 2 ) (2ir) D 5^ D \P - q x - q 2 ). 

We again integrate over q\ and q 3 with the delta functions and introduce 
the on-shell momenta {$ x — 4 3 } and {$ 2 + 4 3 } to obtain 



„{2,0} 
"18 



1 / r + qUj+f^ \ {piMA _ ^} 7a {^ 2+ ^ }7 > 2 > 



MftO V (9s 9 3 

9(A - |p(g 3 ,Pi - 9 3 )|) 9(A - |p(g 3 ,p 2 )|) 
2cj(pi - <f 3 ) p(q 3 ,pi - q 3 ) 2uj(p 2 + q 3 ) p(q 3 ,p 2 ) 



(71) 



Looking at the denominator p(q 3 ,pi — q 3 ) p(q3,p 2 ) it appears that there are 
collinear singularities for g 3 ||p 2 , q 3 \\pi and a soft singularity q 3 — > 0. How- 
ever, the denominator 

(_ r + gg^ + gpn (pi | T { y ^} 7a{ ^ + 4} % | P2 ) (72) 

V (93 9 3 ) / 

vanishes in the collinear regions g 3 ||p 2 and ^Hpi- Thus, this diagram has 
only a soft singularity. 
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Figure 7: Cut diagram for the 2-particle cut diagram with one-gluon 
exchange in the asymptotic region. 



We delegate the explicit evaluation of ajg °^ to the Appendix. Multiplying 
by two to account for both one-gluon exchange diagrams we have the final 
result 

2a\ 2 8 0} = 2jdq 1 dq 2 g 2 Wil\q pl ,q p2 -,q qU q q2 )xA^\q ql ,q q2 ^(P)) (73) 



Pi ~ P2 j 



x / dq 3 Q(A-\p(q 3 ,p 1 -q 3 ,)\)e(A-\p(q 3 ,p2)\)f 2 {pi,P2,q3) 



where again we have not performed the finite f 2 integral analytically and 

<fc(A) = 21o g 2-21og(V). (74) 



2 



The explicit form of f 2 is given in eq. (jl(J2J) . 



4.5.5 3 Particle Cut Diagram 

Let us now turn to the forth part of eq. (|55|). For this term we need the 
dressing factor W^(q p x, q p2 ; q q i, q q2 , g q3 ). Again we use the diagrammatic 
rules of Section I4~4l 
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There are two possible diagrams as the gluon can be absorbed either by 
the quark or antiquark line. The two diagrams are obtained by taking either 
the vertex V\ or V4 and they are symmetric under exchange of momenta. So 
we need only calculate one of them. For the diagram shown in Figure El we 

get 

a\ 1A} = J dq 1 dq 2 dq 3 g 2 W[ 1 \q pl ,q p2 ;q ql ,q q2 ,g q3 ) x A {1) (q ql ,q q2 , g q3 ;>y(P)) 

(_ gPV + <&%±&€ \ e(A-Kg 1 )+^(g 3 )-^(p 1 )|) 
V (?3 9 3 ) / 2uj (Qi) +^(93) -w(pi)) 

\ 2(gig 3 ) 2(g 2 g 3 ) / 

S^iqt + q 3 ~ pi) SP-Qfa - g 2 )(27r) D *W(P q 2 ). 



After integration over gi and g*2 using the delta functions we observe that 
there are collinear singularities q^pi and soft singularities q% — > 0. There are, 
however, no collinear singularities ^3 ||P2- This is expected since the amplitude 
(q q i, q q2 , g q 3', l{P)) has only an integrable square-root singularity for g 3 \\q 2 
and the dressing factor wP(g p i, q p2 ; q q i, q q i, g q z) is regular for g 3 ||g 2 . 

As for the other diagrams we have to multiply by two to take into account 



32 



both pairs of diagrams and we get the final result 

2a\ 1,1} = 2 J dqidq 2 dq 3 g 2 w{ 1 \q pl , q p2 ; q q i,q<p,g<p) (76) 

x A {1) (q qU q q2 ,g q3 ;-f(P)) 
= C F (£J ((| + ; + ^3(A)-c fi ) A(% pl ,q p2 ;j(P)) 

+ (-ie)<% (p 1 |7lp 2 )(27r)( D - 1 )^- 1 )(P-p 1 

x J dq 3 6(A - |p(g 3 ,Pi - &)|)/3(Pi,P2 



9s J 



where 



0s(A) 



7 + 



- 2 log 2 



A 



7t^ 
6 



+ 



A 



-3 + 2 
A 



A 



A 



log 



A 



41ogMl + — -8Li 2 



2 + A 



(77) 



The function / 3 is given in eq. (jl()4j) and does not produce any infrared 
singularity upon integration over dg3. 



4.5.6 An infrared-finite amplitude 

We have now calculated all terms contributing to the amplitude A({q, q}; 7), 
eq. (JHHj), at next-to-leading order. Using eqs. (jHSl EH EH) and (|75|) to assemble 
the amplitude we get 

A({q pl ,q p2 };i)= (78) 
l+C F (g) L(A)+g 2 (A)+g 3 (A)-A+^)A^(q pl ,q p2 n(P)) 

+H'e)^- (p 1 \ 1 a \p 2 )(2^ D ^5^ D - 1 \P-p l -p 2 ) 
x J dq 3 yfi(pi,p 2 ,q 3 )Q{A-\p{ql,p 1 -ql)\)5{^-u(p 1 )-u{p 2 )) 
+/ 2 (pi,P2,?3)e(A- \p(q3,Pi-q 3 )\)&(A- \p(q 3 ,p 2 )\) 

+/3(pi,P2,gs)@(A - |p(g 3 ,pi -93)!) 
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up to order a s in the coupling. The functions g%, g 2 , g 3 are given in eqs. (|b9| 
171)1 and (|77|) and the functions /1, / 2 and f 3 are given in eqs. (fDTJl I1U2|) and 
()104|) respectively. 

We see that this result is completely free of infrared singularities. We 
are only left with some finite A dependent terms, gi and some finite terms, 
fi which will in general need to be numerically integrated. Even though 
the amplitude A({q p i, g^}; 7) depends on A this dependence will disappear 
when we combine the various amplitudes to calculate physical observables. 

4.6 The amplitude A({q(pi), q(p 2 ) : g(p3)}; 7) 

We are now going to calculate the amplitude A({q p i, q P 2, ^3}; 7) given in 
eq. (|56|) . There are only two terms to calculate for this amplitude and there 
is no integration over the final state gluon as it is now a real final state 
particle. 




Figure 9: Cut diagram for the 3-particle asymptotic region with a 
2-particle intermediate state. 

Again we calculate these terms using the diagrammatic rules from Sec- 
tion ^31 Let us start with the diagrams where the gluon is emitted in the 
dressing factor. Figure El shows one of the two possible diagrams, the other is 
exactly the same but with all momenta interchanged. So for both diagrams 
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we have 

dgidg 2 9 M^iVpi, %2, 9 P 3] q q i, q&) x A (0) {q ql , q q2 ; j(P)) 
= (-ie)gTt j (2-K) D 5^- 1 \P-p 1 -p 2 -p 3 ){p 1 \ (79) 

-0(A - \n\)5(^ - a; (pi) - w(p 2 ) - uj(p 3 ) + n) 



& 3 #i+/*3>7 c 
2o;(pi + p 3 )r 1 



2^(p 2 +P3)^2 / 

where we used the notation 

n = p(p 1 ,p 3 )=uj(p 1 )+uj(p 3 )-uj(p 1 +p 3 ) (80) 

T 2 = P(?2,P3) =u{j>2) +^(P3) ~u(p 2 +p 3 ), 

with p defined in eq. (|57|) . 




Figure 10: Cut diagram for 3-particle asymptotic region with a 3- 
particle intermediate state. 



The second contribution is just the usual A^'(q p i, q p2 , g P 3', l{P)) ampli- 
tude. The three external particles of this amplitude do not interact in the 
asymptotic region and so we simply have the diagram as shown in Figure IT0| 
this gives 
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) Qp2, 9p3] Qql, Qq2, QqS) X 

I «^«<pil^ 2(wps) 2( ^ 3) ^ 

(27r) D 5( D )(P- J 9 1 -p 2 -j9 3 ). (81) 
We now assemble eq. (J5T)|) to find 

-4({<?pi,9p2,# P 3};7) = He)#7^(pi|^ 



0(A- InlWy/s - Ufa) -Ufa) -u(p 3 ) + n) 

<y(Vi - w(pi) - a;(p 2 ) - 



^3(^1+^3)7 



2(piPs) 

+ 7 r { ^ + ,^j/ P3 e(A - \r 2 \)5(V~s - - u(p 2 ) - u(p 3 ) + r 2 ) 



2a;(p 2 + p 3 )r 2 

5(y/s - a;(pi) - a;(p 2 ) - a;(p 3 )) |p 2 ) 



7^2+^3) iv. 



2(p 2 Ps) 

(27r)^ D - 1 )(P-p 1 -p 2 -p 3 ). (82) 

This amplitude splits up into two pairs. The first (last) two terms are due 
to the gluon being emitted from the leg p\ (p 2 ). Looking at the first two 
terms shows that for r\ > A the contribution from the asymptotic region 
disappears. We are then left with the normal amplitude, A(q p i, q p2 , g P 3', 7). 
For t"i < A the term from the asymptotic region does contribute and will 
cancel any potential infrared singularities. We can see this by taking the 
limit A — > 0, we have 

^(pi+P3)n -> (pm), 

5(-v/s-^(pi) -w(p2) -w(pa) +n) -> 5(Vs-w(pi) -a>(p 2 ) -o;(p 3 )). 

With these we can see that the terms from the asymptotic region approach 
those of the normal amplitude in the soft and collinear limits, but with the 
opposite sign. So the two terms will cancel in the A — ► limit, leaving 
us with an amplitude that is infrared finite when integrated over the phase 
space. 
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4.7 Calculation of the total cross section 



In the previous sections we computed the two infrared finite amplitudes that 
contribute to the process 7*(P) — > 2 jets at next-to- leading order. In this 
section we would like to check our results by computing the total cross section, 
starting from the infrared finite amplitudes, eqs. (J78j) and (JB2J- Of course, 
we have to recover the well known result, eq. (|54|). 

Let us stress that the idea of our approach is to compute the amplitudes 
numerically and perform the phase-space integration also numerically. It is 
for the sole purpose of checking our results and facilitating the comparison 
with eq. (J53)l that in this section we compute the total cross section analyti- 
cally. 

Usually a non zero value of A would be chosen for a numerical calculation 
and we would expect all A dependence to cancel between the contributions 
of the two amplitudes (squared) to the cross section. Here though to simplify 
the analytical calculation we will take the limit A — > 0. In this limit even 
the infrared-finite amplitudes are proportional to a four- dimensional delta 
function and we can use the standard procedure to obtain the total cross 
section from the amplitudes. However, since the amplitudes are singular for 
A — ► we must be careful in taking this limit and leave it until the end of 
the calculation. The f n finite terms of eq. (J78j) which we would usually have 
to calculate numerically will all go to zero in this limit. This simplification 
occurs because the region of integration shrinks to zero as A — > and as 
these terms are finite they can no longer give a contribution. 

We now use our infrared finite amplitudes eq. (|55|l and eq. (J55J) instead 
of eq. (J45|) and eq. (J46|) and square them in the usual way to obtain 



Here we integrate eq. over the two particle phase space and eq. (JHo"|) over 
the three particle phase space. 



(83) 



where 




(84) 



(85) 
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First we rewrite the three-particle final state amplitude, eq. (|82|) . in a 
more convenient form 

<A({q p uq P 2,g P 3hj) = (-ie)gT t a j (pj|^ 



<J(£? + n) + ^vri ^ Q(A _ | ri)) 



2o;(p*i+p 3 )ri 2(pip 3 

«J(E)e(|n|-A) 



^1+^)7° 



2(piPs) 

\ 2uj(p 2 +P3)r 2 2(p 2 P3j / 



7^2+^) ^3 



5(£)6(|r 2 |-A) |p 2 > 



2(P2Ps) 

(2 7 r) D 5^ l \P-p 1 -p 2 -p 3 ), (86) 

where .E = \fs — to{p\) — ^(p 2 ) — uj(p 3 ). Taking eq. (JHHj) we then square it in 
the usual way and sum over the gluon polarizations using 

2. P3 ps - 9 + te) 

where p 3 = (o;(p 3 ), — p 3 ). This is because the amplitude is no longer gauge 
invariant as we are using dressed states. At this point we drop any terms 
multiplied by 0(A — \ri\) or @(A — |r 2 |). These terms are finite and therefore 
can be shown to go to zero in the A — ► limit after we have performed the 
three particle phase space integral in a similar way to the f n terms. 

After integrating one of the phase space integrals using the delta function 
we are left with 

\A({q pl ,q P 2,9 P shl)\ 2 = 4C F (^?- 2D J (87) 
" "° , 2y 23 - y 13 (y l3 + y 23 f 



f 2 dy 13 f dy 

JO Jl— I/1S 



/0 J i-yia y 13 [yi 3 + y 2 3) 

j ft A Z/23 + I/13Z/23 + 2yi 3 (y 23 - 1) 

"' 23 ~t i \2~ 



/ 2 ^2/13 / 2 ^2/2u 

Jo Jl -yw y 2 3 [yi3 + 1/23 

f 1 -^ , ff f 2-y 23 2 -yi3 



i-tna V 2/13 1/23 (yi3 + ?/ 23 )' 
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where we defined 

* - (88) 

We perform the final two integrals and then prematurely take the A — > 
limit everywhere except in the log(A) terms, as these diverge in this limit. 
The log(A) terms will cancel later in the final result. This then leaves 



-<«=Mv)(^ 4 +Ht) +bg2 (t)) (89) 



1 „ /a„\ (h , . 3 , /A\ , 9 /A N 

where is the total Born cross section as given in eq. 

Now we calculate |^4({g p i, q P 2}] Again we take the A — > limit early 
except for the log(A) pieces of the g n terms in the finite part of eq. (J7HJ). As 
stated before the f n (pi,P2, 93) terms go to zero and so we have, 

l^({g P i,g p2 };7)| 2 = \^ \q pl ,q p2 ;i(P))\ 2 (90) 

After integrating over the two particle phase space we get 

^ - ( l + C F * + log4 - I log ($) - log 2 ($))+ 0(al 



(91) 



a \ 7T V 2 2 

Putting eq. (|8*9*|l and eq. (|9T|) together gives finally 

^ = ( 1 + (T)i c ' + ^)- (92) 

We have recovered the well known result for the total 7 — > qq cross section 
and all the A dependence of the amplitudes has disappeared including the 
log(A) terms, justifying our taking of the A — >• limit early. 



5 Summary and outlook 

We have presented a method on how to construct infrared finite amplitudes 
and applied it to the case of e + e~ — ► 2 jets at next-to-leading order in the 
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strong coupling. The idea is to separate from the Hamiltonian a part that 
describes the asymptotic dynamics. This asymptotic Hamiltonian is then 
used to asymptotically evolve the usual states of the Fock space. In this 
way we construct dressed states, eqs. and (j23), such that the transition 
amplitudes between these states are free from infrared singularities. 

Contrary to most of the previous work done in this field we are not so 
much interested in obtaining all-order resumed results taking into account 
soft emission of an arbitrary number of gauge bosons from external partons. 
Our aim is to construct dressed states explicitly order-by order in perturba- 
tion theory and use them to do explicit calculations. In this paper we have 
done this for a particularly simple final state up to next-to-leading order. In 
the future we would like to expand this to more complicated external states 
and higher orders. 

The reason that we cannot obtain all-order results is that we include 
the collinear singularities as well. In non-abelian theories these singulari- 
ties cannot be avoided. The additional complications due to the collinear 
singularities make it impossible to obtain exact solutions to the asymptotic 
dynamics. Collinear singularities have been considered previously 
but to the best of our knowledge the amplitudes presented in this work are 
the first infrared-finite amplitudes for a realistic scattering process in QCD. 

As for the standard approach, physical cross sections obtain in general 
contributions from more than one partonic process. However, in our case 
all these contributions are separately finite. They depend on a parameter, 
A, that determines the precise split of the Hamiltonian into an asymptotic 
Hamiltonian and the remainder. The result for any physical quantity is 
independent of this parameter as long as it is smaller than any experimental 
resolution. For any finite value of A the amplitude contains a part that is 
not proportional to an energy conserving delta function which represents the 
spread of the initial wave packet due to the asymptotic evolution. 

For any physical cross section at any order in perturbation theory we will 
get the same answer using the standard cross-section method or infrared- 
finite amplitudes. Thus, one might wonder what has been gained using this 
approach. Apart from the conceptional benefit that the ^-matrix between 
dressed states is well defined there are also practical advantages. First of 
all, the avoidance of infrared singularities facilitates the use of numerical 
methods. This might not be apparent in the approach we have taken. In 
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fact, using eq. (|28)l to split the infrared finite amplitudes into separately 
divergent factors still requires us to use an infrared regulator (dimensional 
regularization in our case) and revert to analytical calculations. However, 
since the final amplitude is infrared finite it is feasible to compute it directly in 
a numerical way, avoiding the split into separately divergent pieces. Once the 
amplitudes have been obtained, the integration over the phase space is trivial 
and no sophisticated method is needed. This also opens up the possibility of 
combining fixed-order calculations directly with a parton shower approach. 

Needless to say that the explicit example we considered, e + e~ — > 2 jets 
has many simplifying features. To start with, the non-abelian nature of QCD 
does not really enter. Secondly, we only considered the amplitudes at next- 
to-leading order. Furthermore, the initial state does not interact strongly. 

The last point simply results in the fact that there is no need to dress 
the initial state. While this is a simplification concerning the amount of 
computations to be performed, there is no conceptual problem associated 
with more complicated initial states. If the initial state contains hadrons a 
physical cross section is obtained by folding the partonic cross section with 
parton densities. In the conventional approach these parton densities are 
associated with the probability of finding a certain partonic state within a 
hadron. In our case, we would have to use modified parton densities that are 
related to the probability of finding a certain dressed state within a hadron. 
Thus the global analyzes of extracting the parton densities would have to be 
modified and repeated. 

The fact that the non-abelian nature of QCD does not really show up in 
the explicit example we considered results in a particularly simple asymp- 
totic Hamiltonian. In fact, the asymptotic Hamiltonian we use involves only 
quark-gluon interactions and is basically the same that was used many times 
previously ^3]. Again, this results in a technical simplification of the com- 
putation and facilitates the explicit construction of the asymptotic Hamil- 
tonian. In more complicated examples the full non-abelian structure of the 
asymptotic Hamiltonian will enter the problem and its construction will be 
much more involved. However, the only crucial feature is that the asymp- 
totic Hamiltonian reproduces the full asymptotic dynamics, i.e. it has to 
reproduce the soft and collinear behavior of the full theory. There are no 
further requirements and the construction of dressed states presented in this 
paper can be taken over directly. However, it is clear that the construction 
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used so far is rather cumbersome. In order to exploit the advantage of the in- 
frared fmiteness a systematic numerical approach should be developed. This 
will become particularly important if the method is to be extended beyond 
next-to-leading order. 



Appendix 

In this Appendix we give some details concerning the evaluation of the dia- 
grams mentioned in Section 1431 

We consider first with the self-interaction term a\ 2 ' °^ of Sect ion 14.5.31 We 
start from eq. substitute^ — / for {j^ — j/ 3 }, where r = (r , 0) with 
r o = p(Q3,Pi — Qa) and then expand the numerator to obtain 



{2,0} He) g 2 



15 



T^(27T) D 6^(P- Pl -p 2 ) Jdq 3 e(A-\r \) (93) 



4 



(Pl|7"|P2) 



w(pi)w(^i - q 3 )r% 



This expression contains infrared singularities coming from the region where 
#3 is soft and/or collinear to p\. In order to evaluate the expression, eq. (jHSJ) 
we choose to parameterize the momenta in the center-of-mass frame. The 
momenta are all on-shell and are defined as 

P = Vi(l,0,0), (94) 

Pi = %(i,o,i), 



P2 



%(i,o,-i), 



q 3 = |lz(l,yr^7e T ,y), 



where is the null vector in a (2— 2e)-dimensional space, is a unit vector in 
the (2— 2e)-dimensional transverse momentum space and we have < z < 00, 
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— 1 < 2/ < 1. The singular limits are then given by the limits z — > for soft 
singularities, y — > 1 for qs\\pi singularities and y — > — 1 for qz\\p2 singularities. 
As the asymptotic region does not conserve energy we find that the upper 
limit of z goes to oo. This would suggest the possibility of UV singularities in 
the asymptotic regions. However we will see that the function will restrict 
this upper limit to a finite value, removing the need to renormalize these 
regions. The integral measure is given by 

(95) 




s ) 2(27r)3 - 2 e f z^(l-yr(l + yrdydzdn {2 _ 2e) 

where we have three separate integration regions for the z and y integrals, 

2 (2 + A) , 2 -A 2 

< z < — J — with -1< y < , 

- - A(l-y + A) ~ y - 2' 

2- A 2 

< z < 1 with — - — < y < 1, 

2 (2 + A) , 2 -A 2 

1< z< — — i J — with < y < 1. 

" -A(l-j/ + A) 2 ~ y - 

The infrared singularities are in the first two regions whereas the last region 
will give a finite contribution. The remaining angular integral is given by 

<aV &) = f ^y- (96) 

We now turn back to eq. (}93|) and notice that the infrared singularities qs 
soft and/or collinear to pi come from the region z = and y — 1 but not 
y = —1. We use the subtraction method to isolate these singularities and 
evaluate them analytically. Writing the integrand schematically as a function 
F(z, y) we write 

F{z,y) = (F(0,y) + F(z,l)-F(0,l)) (97) 
+ {F(z,y)-F(0,y)-F(z,l)+F(0,l)). 

The first term contains all the divergent pieces whereas the second term will 
give a finite contribution upon integration over dg3. Applying this method 



43 



to eq. (|93|) we obtain 

= H^ 2 ^-(^] n/n _V, f ~ Pi ~ fi) (98) 




(Pi|7>2> / ^( 2 -2e) dydzz 1 ^ (1 - y)- (1 + ?/) 



f(2-D) (2z-l-y) , , 



where 



/l(Pl S P2 S ?3) = (99) 
w(Pl) 2 A (PlftO (Pl9s)(P2«3) / 



1 - 



2a; (pi - g 3 ) r \ a; (pi) r w(pi)r (g 3 g 3 ) 
^(pi) 2 / w(g 3 ) (p 2 g 3 ) , \ 




2(pig 3 ) I a; (pi) w(g 3 ; 



2 



(100) 



Integrating the singular terms and expanding around e = we obtain eq. 
Note that the D in the first term of eq. (|9"8j) arises from the 7-matrix algebra. 
Thus we write it as D = 4— 2e+Cfi2e to obtain the expressions in conventional 
dimensional regularization (cr = 0) and in dimensional reduction (cr =1). 

Let us now turn to the evaluation of a\l' ^ needed in Section 14.5.41 We 
start with the expression eq. (J7TJ) and proceed in the same way as for a| 2 ' ^. 
We introduce the on-shell momenta {p\— ^ 3 } =i>\— $3+ f and {$ 2 + ^ 3 } = 
^2+ #3- f, where r' = (r' ,0) with r' = p(g 3 ,p 2 ) = u(q 3 ) +u(p 2 ) -u(p 2 + q 3 ). 
In order to proceed we subtract the soft singularity in eq. (fTTj) and add it 
back to produce an integrand that results in a non-singular term. In the soft 
limit the D- dimensional delta function becomes the usual 8^ D '(P — pi — p 2 ) 
which can be pulled out from the integral. 

We use the same momentum parametrization as for the self-interacting 
case, but because of the extra 6 function the integration ranges change to 

2 (2 + A) , A 

< Z < — t± —rr With -1< V < — , 

~ ~ A(l -y + A) - y ~2' 

2 (2 -A) , A 

< z < -— '— with — < y < 1 

" " A (1 + j/ -A) 2 ~ y - 

The remaining angular integral is as given in eq. (J96|) . 
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Using this momentum parametrization and expanding around the soft 
region gives 



,{2,0} 



2 N 

s) v 2(2vr) 

l-2e 



(Pi\l a \P2) I dnp.. 2e) dydzz L -"(l-y)- e (l + y)- t 

(^(v^ - w (pi) - ^)) + /2(pi>P2, ga)) (101) 



with 



f2(pi,P2, ?s) = TTT^ Jt^/^ i ^\ / ( (P1P2) + (Pi9s) - (p^a) 

(pm)(p2g3) 1 / {pigs) 3 L 
2u(p 1 y 2(^3) Wn) 2 \ 

(P2qs) 3 (, r' \ + / ( r 



r 


) 


2w(ft) 






r r£, 




u;(pi) 2 


- w(p 2 ) 


- r + ^d) 



- |^«(V5-w(ft)-w(fi)). (102) 

Upon performing the integration of the singular terms explicitly and expand- 
ing in e we get eq. (|73J) . In this case the expression is the same in conventional 
dimensional regularization and dimensional reduction. 

Finally we turn to the evaluation of a} 1 ' 1 ^ needed in Section I4.5.5| pro- 
ceeding as in the previous cases. We subtract the soft and collinear singular 
parts and integrate them analytically. In both limits the D-dimensional delta 
function takes its usual form 5^ D \P — p\ —p2)- Thus, the delta function is in- 
dependent of the integration variables and can be taken outside the integral, 
as in the one-gluon exchange terms. 

We can use the same momentum parametrization and integration regions 
as the self-interacting case as we have the same G function in both cases. 
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Taking the 2 — > and y — > 1 limits of the above terms we obtain 



= H^^(^) ^^(2 7 r)^- 1 )(P-p 1 -p 2 ) 

(Pi|7>2> / dn^ 2e) dydzz l - 2e {l-y)- e {l+y)-* (103) 

/4-4z+(D-2)z2 A 
x I 2 ^ _ ^ o(vs-w(pi) -w(p2)) + / 3 (pi,p2,fe) I 

where 

_ (p,^) f 2 + _r° 1 + (M%gj f 2 + Jf±_ + '0 ^ 



w(pi)y (g 3 g 3 ) V ^(pi) MpO 



2 



+ ^(Pi) 2 f 2(p!p 2 ) + (pig 3 ) f 2 + r ° - ( ff2g3 ) ^ 

-2^3) +2^0-^(2+ r ° 



(9393) V w (p0 

+ ^( 1+ 2^))))^-( P -)-^)- r ») 

1 UPi^sX^gs) w(pi)(Pi93) 2(p 1 g 3 )w(g 3 ) 2 
(Pi?a) 2 



2(P2<?3V(<?3) 2 (pigs) 



|<f(v^ (104) 



Integrating the singular terms with D = 4 — 2e + c^2e and expanding around 
e = we obtain eq. (|75|l . 
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